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Abstract 

In this work we investigate the role of the symmetry of the Lagrangian on the existence 
of defects in systems of coupled scalar fields. We focus attention mainly on solutions 
where defects may nest defects. When space is non-compact we find topological BPS and 
non-BPS solutions that present internal structure. When space is compact the solutions 
are nontopological sphalerons, which may be nested inside the topological defects. We 
address the question of classical stability of these topological and nontopological solutions 
and investigate how the thermal corrections may modify the classical scenario. 



PACS number(s): ll.10.Lm, 11.27.+d, 98.80.Cq 



1 This work is supported in part by funds provided by the U.S. Department of Energy (D.O.E.) under 
cooperative research agreement #DF-FC02-94ER40818. 



1 



1 Introduction 



In recent years much investigation on topological defects have been done with applications 
on diverse fields as for instance cosmology and particle physics [|], 0. In the present 
paper we are interested in systems of coupled scalar fields where solutions like defects 
inside defects may appear. Such situation was first investigated in the work of Witten 
||, within the context of superconducting strings - see also the work of Lazarides and 
Shafi [Q and of MacKenzie ||. Several aspects of the possibility of topological defects 
being nested inside topological defects are considered in Refs. || [7], §, |9|, pTQ|| , and in 
Refs. JO], [12], |13|, [13], [15|] other related issues are also investigated. 



To introduce defects inside defects one needs to consider systems of at least two fields, 
one of the fields being able to describe a specific defect, which will respond for nesting the 
other defect in its interior. The system of two fields should then be able to accommodate 
some defect and, in the core of this defect, be reduced to a system of the second field that 
is still able to generate the other defect. Since these defects are usually topological, and 
since topological defects appear in systems presenting spontaneous symmetry breaking, 
the complete system usually presents Gi x G 2 symmetry, for instance U(l) x U(l) in the 
case of the superconducting strings examined by Witten or Z 2 x Z 2 domain walls, as 
considered for instance in the recent works |7|, |], [TDfl . 

Although systems of two coupled fields engendering G\ x G 2 symmetry are known to 
support defects inside defects, we may wonder if defects with internal structures appear 
after changing the symmetry of the system. Stated differently, we may ask how the 
picture of defects inside defects modifies when one changes the symmetry of the underlying 
Lagrangian, enlarging or reducing the G\ x G 2 group. This issue is of general interest and 
may find direct applications in cosmology and in condensed matter. In the present work 
we address this and other related questions in systems of two real scalar fields, described 
via the Lagrangian density 

c = l -d a <pd a <p + l -d aX d a x - U(<f>, x) ■ (i) 

Although we are working in 3 + 1 space-time dimensions we can search for static solutions 
in the form (p = (f)(x) and \ = x( x )i depending on only one (x 1 = x) of the three spatial 
coordinates. In this case the equations of motion are given by 



2 /u ' ( 2 ) 



fj> dU 
dx 

fx = dU 

dx 2 dx 
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We are interested in a special class of systems, in which the potential is determined by 




(4) 

It is defined in terms of a smooth but otherwise arbitrary function H = H(<f), x) and now 
the equations of motion become 



d 2 



dx 2 

fx 

dx 2 



H^Hjjcf, + H x H x< f, , (5) 
H^Hfa + H X H XX , (6) 



where stands for dH/d<f), and so forth. The advantage of considering such potential 
is that the energy of static configurations is minimized to the value 

E B = \H[<f>(+oo),x(+oo))-H[<f>(-oo),x(-oo)}\ , (7) 

which depends just on the asymptotic values of the field configurations. Furthermore, the 
second-order differential equations of motion are now solved by field configurations that 
obey the first-order equations 

This result can be shown [|TTJ to follow the idea of Bogomol'nyi, Prasad and Sommerfield 
(BPS) |n| in the case of coupled real scalar fields. This means that the potential in Eq. (f|) 
can be considered as describing the bosonic sector of a larger supersymmetric theory, so 
that supersymmetric extensions of the Lagrangian ([!]) can be readly constructed, as was 
considered for instance in Ref. [ID, [T^J as well as in several other investigations [JT^J, FJ, ITS]. 



The present work is organized as follows. In the next Sec. |2] we search for defects 
inside defects starting with a general potential. There we show that defects may engender 
internal structure only when the system presents parity symmetry of the type Z2 x Z2. 
Once this is stablished we write down explicit solutions, examining both the second- 
order equations of motion and the first-order equations that give BPS solutions. We 
also consider the possibility of dealing with compact space, to obtain periodic solutions 
similar to the sphaleron solutions discussed by Manton and Samols [H]]. Such solutions 
are known to be useful for discussing violation of the barionic number, and may contribute 
to applications in cosmology. In Sec. |3] we investigate classical or linear stability of the 
solutions introduced in Sec. 0. Although some of these solutions were already discussed 
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in Refs. || ||, a complete investigation of their stability is still missing and will be done 
in this work. In Sec. |] we introduce finite temperature contributions to investigate how 
the thermal effects may change the classical scenario. We end this paper in Sec. |5|, where 
we present comments and conclusions. 



2 Symmetry and defects inside defects 

We start with the simplest case in which a theory presents a Z 2 parity symmetry with 
nontrivial vacuum states for a single field <fi. For 

U(4>) = ±H}, (10) 

the standard </> 4 model is obtained with the function 

H(<p) = X(^-a 2 <p). (11) 

A natural generalization for two coupled fields, in a renormalizable theory in (3+1) di- 
mensions can be written as 

H(<f>, x) = A(^0 3 - a 2 0) + W>X 2 + W> 2 X + <\x" - b 2 X ), (12) 
where A, /i, v, a, a and b are real parameters. Here we have the general potential 
U(<f>, X ) = ^(A 2 + z/ 2 )0 4 + 2z/(A + /i)0 3 x + (A/i + 2z/ 2 + az/ + 2 yU 2 )0V 

+ 2fi(u + a)<p X 3 + \(^ + * V - + °b 2 v)<p 2 - (Aa 2 /i + a 2 b 2 ) X 2 

- 2(Aa 2 z/ + a& 2 /i)0x + ^(AV + a 2 6 4 ) . (13) 

We now impose that the theory presents Z2 x Z2 symmetry. To ensure invariance of the 
Lagrangian under the independent transformations — > —(f) and x ~^ ~X we se t 

u{\ + n)=0, (14) 
M (i/ + a)=0, (15) 
Az/a 2 + (T/i6 2 = 0. (16) 

There are four possibilities of satisfying the above conditions; they are: 

(i) v = a = 0, (17) 

(ii) A = /i = 0, (18) 

(m) 1/ = /x = 0, (19) 

(iv) A = -w v = -ex; a 2 = -b 2 . (20) 
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The first possibility gives a model for which the function H(<fi, X ) reads 



X) = - a 2 <P) + fi<Px 2 , (21) 

which was discussed recently in the literature J7|, ||. This is a tipical example, shown to 
present topological solutions of the type of defects inside defects. The second possibility 
leads to a model that is described by 

H 2 ( ( p,x) = ^X 3 -b 2 x) + u ( p 2 x, (22) 

which is equivalent to the previous one, as can be seen by simply exchanging the fields 
and x an d redefining the parameters a, v and b. The last two possibilities give uninterest- 
ing models, at least from the point of view of defects inside defects: the third possibility 
leads to a system of two decoupled scalar fields and the fourth one severely restricts the 
number of nontrivial vacua, forbidding the presence of defects inside defects. 

2.1 Topological solutions 

We consider the equations of motion (H) and (|6]) for the general model described by 



H((f),x), given by eq. (p"2"D . We have 
d 2 6 



dx 2 

fx 

dx 2 



[A(0 2 - a 2 ) + /i X 2 + 2W> X ](2A0 + 2v X ) 
+ [2^0 X + W> 2 + °{X 2 - b 2 )](2fi X + 2W>) , (23) 
[A(0 2 - a 2 ) + fix 2 + 2u<Px](2fx X + 2v<j>) 

+[2fuf> X + "<f> 2 + <X 2 ~ b 2 )}(2^ + 2ox) ■ (24) 



To find defects inside defects we have to search for pairs of solutions like (0,0) and (0,x)- 
For the first pair we set X — an d then the above equations reduce to 

= [2A 2 (0 2 - a 2 ) + 2z/ 2 2 - 2avb 2 \<\> , (25) 

dx 2 

and 

v{\ + n)(p 2 - (Az/a 2 + fxab 2 ) = . (26) 
For the second pair we have = and then Eqs. (|23"D and (|24]) read 

g = [2a\ X 2 - b 2 ) - 2^Xa 2 + 2^\ 2 ]x , (27) 

5 



and 

Hiy + a)x 2 - (Az/a 2 + fiab 2 ) = . (2? 
Eqs. (PB|) and (123) show that solutions like defects inside defects are only possible for 



u(\ + fj) = 0, + or) = 0, Xua 2 + fiab 2 = 0, (29) 

Interestingly, these conditions are exactly the conditions (|14[)-(|T6D , needed to ensure the 
Z2 x Z2 parity symmetry of the system. This means that only the generating functions 
Hi(cj),x) and H 2 (<f),x) ma y admit solutions describing defects inside defects. 

The explicit solutions can be found by solving the remaining equations. For the system 
described by Hi we take v = o = and rewrite Eqs. fl25|) and (p7|) as: x — and 



' h " 2A 2 (0 2 -a 2 )0, (30) 



dx 2 
or = and 

g = 2/ , 2 ( X 2 -ra 2 ) X , (31) 
where we have set A//i = r. There are pairs of solutions 

(f)(x) = atanh(Aax) ; x — , (32) 

and 

x(x) = a\pr tanh(/^ Vraa;) ; = 0. (33) 

These solutions present the feature of vanishing at its own core, that is, for x —>■ we 
get (f>(x) — > and x( x ) ~^ 0- This is important for introducing defects inside defects in 
3 + 1 dimensions, since in this case we can choose = <p(x) and x = x{y)i sa Y) to make 
topological defect appear in the core of topological defect. The issue of choosing the field 
to host the second field depends on the parameters of the system, and follows as in the 
former investigation || - see also Ref. [0 for similar considerations in 2 + 1 dimensions. 



Similarly, we can also find solutions like defects inside defects in the other system, 
defined via ^2(0, X)- I n this case the pairs of solutions are given by 

<j,(x) = 6y^tanh(z/y^6x) ; x = , (34) 

and 

x(x) = 6tanh(cr6x) ; = 0. (35) 
Also they may be used to build defects inside defects in the 3 + 1 dimensional system. 
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2.2 BPS solutions 



Let us now consider the first-order equations (§) and @ for the general H given by 
Eq. (|P2|). Here we have 



We set x — to get to 



and 



A(0 2 -a 2 )+/ix 2 + 2z/0 X , (36) 
= a( X 2 - b 2 ) + z/0 2 + • (37) 

|j = A(0 2 -a 2 ), (38) 



V( f) 2 -ab 2 = . (39) 
We choose z/ = o = to get to the pair of solutions 

4>(x) = atanhAax ; x — ■ (40) 

Analogously, in the case <fi = we have 

5r = °(x 2 - h2 ) . ( 41 ) 

and 

/i X 2 - Aa 2 = . (42) 

We impose A = /x = to obtain 

x(x) = fotanhcrfer ; = 0. (43) 

These are BPS defects, since they arise from the first-order equations @ and (§). 

It is interesting to realize that the above BPS solutions cannot be solutions of a single 
system, because the conditions v = a = and A = /i = cannot be implemented 
simultaneously in the general H given by Eq. {[T2]). This conclusion is in agreement with 
results of dynamical systems, since one knows from unicity of solutions that any two orbits 
never cross each other in configuration space. And we recall that the pair of first-order 
equations (|36D and ([37]) can be seen as a dynamical system. However, we can have the 
BPS defect (|40| ) in the system defined via Hi((f>, x), and the BPS defect ( f43D in the other 
system, defined via H 2 (4>,x)- 

The above solutions have been investigated before || |9|] , but here we have a stronger 
result, showing that systems engendering solutions like defects inside defects must present 
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^2 x Z2 symmetry. Furthermore, only one of the solutions can be of the BPS type. Since 



the BPS solution is stable [IT], we still have to investigate stability of the other solution 
to talk about stability of defects with internal structure. 

For completeness let us investigate the presence of other BPS solutions in the system 
defined by H±((f>, x), for instance. In this case the first-order equations are 

d ± = A(0 2 - a 2 ) + ^ , (44) 
ax 

5r = 2 ^X- (45) 
ax 



An interesting pair of BPS solutions was already obtained in [11|. It is 



4>(x) = — a tanh(2 fiax) , (46) 
x(x) = ±ay/r — 2 sech(2yuas) . (47) 

It is valid for r = A//i > 2 and obeys 

t 2 + ^L=X 2 = a\ (48) 
\/r - 2 

which defines a semi-ellipsis in configuration space. We recall that at r = 3 the amplitude 
of the two above solutions degenerate to a single value, the orbit changes to a semi-circle 
and the corresponding stability can be implemented analytically J33 . 



This system presents peculiar behavior at two other values of r. For r = 1 the 
symmetry changes from Z 2 x Z 2 to Z 4 , but now the system degenerates into two uncoupled 
systems of one field H. For r = — 1 the vacuum states at the x ax i s desapear, and 
although there is still the pair of BPS solutions given by Eq. (|40|), it is possible to show 
that solutions like the above one, describing a finite orbit connecting the two vacuum 
states (—a, 0) and (a, 0) with non-vanishing x cannot be present anymore. The proof 



follows after recognizing that for r = — 1 the first-order equations (H) and ([EJ) implies 
that 

X (f - \ X 2 ~ a 2 ) = . (49) 
This restriction shows that x = or 

<P 2 = a 2 + \ X 2 (50) 

The first condition x = gives exactly the case of the straight line connecting the vacua 
(— a, 0) and (a, 0). The other condition Eq. (p0|) shows that for x no finite orbit can 
connect the vacuun states anymore. 



S 



2.3 Periodic solutions 



Let us investigate the same systems given by Hi and H 2 , but now searching for periodic 
solutions of the second-order equations (||) and (|3|) that appear when the line becomes 
a circle, a compact space. In this case the equations of motion do not change but the 
boundary conditions are now periodic ones, to respect the topology of the circle. The 
solutions here are like the sphalerons solutions introduced in |L7|| in the case of a single 



field - see also ||TE|, [HI and references therein for further details. For generality we shall 



assume that the fields and x are periodic according to 

<f)(x) = <f)(x + Li) , (51) 
X (x) = x(x + L 2 ) . (52) 

Evidently, when these fields are immersed in a one-dimensional space we should impose 
that Li = L 2 , but in higher dimensions we can search for solutions with = 0(x) and 
X = x{y), f° r instance, and so we can allow Li ^ L 2 and this is the case we shall consider 
below. Here the solutions will be given in terms of the Jacobi elliptic function sn(x), 
which has period 4K(ki) determined by the elliptic quarter period K(ki) which obeys 

(fce(0,i),i = l,2) @ 

ivr = K(0) < K(ki) < K(l) = oo . (53) 

For the model Hi we get the solutions: x = and 

<fi{x) = akib(ki) sn (a\b(ki) x) , (54) 

and also = and 

x(x) = ak 2 y/rb(k 2 ) sn (afi^/r b(k 2 ) , (55) 

where 



The stability of these solutions will be studied in detail in the next section. Analogous 
solutions for H 2 can also be found by just changing — > x and redefinig the parameters 
conveniently. 

In the limit ki — > 1 the radius of the circle goes to infinity, sn(b(ki)x) — > tanh(x), 
and the above solutions recover the corresponding basic kink solutions we have already 
found in the previous subsections. Here we notice that in two or more space dimensions 
we may set ki — *■ 1 or k 2 — *■ 1, to make periodic or topological defects to host topological 
or periodic defects, allowing for new pictures of defects inside defects. As one can check, 
none of the above periodic solutions of the second-order equations of motion solves the 
first-order equations, despite the fact that the limit ki — > 1 recover all the topological 
solutions including the BPS ones. 
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3 Linear stability 



In this section we investigate the stability of the solutions presented in the former Sec. ES. 
Evidently, a general study of the stability of static solutions of nonlinear equations is a 
highly nontrivial matter, but the simplest situation known as classical or linear stability 
may be performed and will be analysed with some detail below. Firstly we recall that the 



solutions (|3^) and (|35| ) also solve the first-order equations, and so are linearly stable |L1 
All other solutions should be investigated with respect to their stability. 
We start by writing the time-dependent fields in the usual form [[H], |22fl 

(f)(x,t) = (j)(x) + rf(x, t) , (57) 
X (x,t) = X (x)+C(x } t) , (58) 

where r](x,t) and COM) are the time- dependent perturbations, small when compared to 
4>(x) and x( x )i the static solutions of the second-order equations of motion. Since we are 
considering fluctuations about static solutions, we can always write 

r)(x,t) = ^2r) n (x) cos(w n t) , (59) 

n 

COM) = XXnOc) cos(tM) • (60) 

n 

Substituting these expressions into the time-dependent Euler-Lagrange equations of mo- 
tion we get the following Schrodinger equation 

+ M]M = «* M , (61) 



dx J y Cn / \ Cn 

where 1 is the 2x2 identity matrix and 

M -(Z £)■ (62) 

The derivatives of the potential U(4>, x) should be calculated at the classical values <j>(x) 
and x(x), and we recall that stability means absence of negative 

The standard analysis of stability usually considers diagonalizing the above matrix 
M, but this may make the problem untractable analytically. We illustrate this situation 
by following the lines of Ref. f23|, where instead of considering the matrix involving 
derivatives of the potential U(<f>,x), one considers another one, involving derivatives of 
H((p, x)- This procedure follows naturally from the fact that the first order equations can 
be used to simplify investigations concerning the second-order equations, at least when the 
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classical configurations also solve the first-order differential equations. For BPS solutions 
this is implemented after introducing the first-order operators 



±l4~ + m 

dx 



where m is the matrix 



m 



\ H x<t> H xx ) 
so that the stability equation fl6"T| ) can be rewritten as 



d 2 dm 2 

dx 2 dx 



Vn 

v Cn . 



S-L- S— 



Vn 
Xn , 



10: 



Vn 
Xn. 



(63) 
(64) 

(65) 
(66) 



In ( |65|) the Schrodinger-like operator presents the form S+S-. However, in Eq. (|63[) the 

-1 



first-order operators are such that S± = S^. Thus we can write 



w. 



< n\S+SXn >=< n\S J _SXn >= \\S-\n > 



> 



where \n > stands for the (orthonormalized) state 



\n > 



Vn 
Xn. 



(67) 



(68) 



This shows explicitly that the Schrodinger-like operator in Eq. (|65 ) is positive semi- 
definite, and so the corresponding eigenvalues are non-negative - see Ref. [[II]] for 
further details. The BPS solutions are stable, but if we want to find the eigenvalues 
explicitly we need to diagonalize the corresponding potential. This diagonalization for 
BPS solutions can be given in terms of the diagonalization of the simpler matrix m. In 
this case the elements of the diagonalized matrix are 



Ah 



- t (H H + H xx ) ± -\[{H 4 



H ) 2 
11 xx) 



AH, 



4>x •> 



and the corresponding Schrodinger-like equations are given by 



d 2 



dx* 



+ 



d\± 
dx 



(A, 



2 ± 

TV 71 
n In 



(69) 



(70) 



The problem of solving the above equations analytically is now related to the presence 
of the square root in the potential of the corresponding Schrodinger equations. To mantain 
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the investigation analytical, we should get rid of the square root in Eq. (|69|), and this can 
be done by imposing one of the following conditions 

(i) : = , (71) 

(ii) : H„ = H XX , (72) 

(iii) : H H H XX = H^ 2 . (73) 

For the general model, the above condition (i) can be satisfied with the following possi- 
bilities 

(ia) : ji = v = , 

(ib) : v — and % = , 

(ic) : fi = and = 0. 

The possibility (ia) must be discarded since it decouples the two fields. The other two 
possibilities (ib) and (ic) can be considered and here <fi in case (ib) and x i n case (ic) 
may describe topological domain defects of the BPS type, but in these cases no defect is 
allowed to appear inside it. 

The second possibility of eliminating the square root from the stability equation is 
= H xx . Since the generalized model presents H^ = 2A0 + 2z/x and H xx = 2fi(p + 2ax 
we see that the condition (ii) is obtained by imposing A = fi and a = v. In this case we 
have 

H s = fi Q</» 3 - a 2 ^ + ^x 2 + H\ + v Qx 3 - b 2 X ) ■ (74) 
This function can be written in the form H s = Hi + H 2 , where 

Hi = juQ^-aV) +W> X 2 , (75) 

H 2 = v Q X 3 - b 2 X ) + vxtf ■ (76) 

In this case, however, we follow || to see that each Hi (i = 1, 2) can be written in terms 
of 4>± = 2 -1 / 2 (0 ± x)-, i n a form that decouple the two fields, and so also does the full 
system described by H s . Thus we see that the general model does not give any system of 
two coupled fields under the additional condition H^ = H xx . 
The third and last condition (iii) implies 

Xficf) 2 + avx 2 + (Act + H0X = ^4> 2 + /^Y + , (77) 

which is satisfied when we identify A/i = v 2 and av = fi 2 , which also implies Act = fiu. 
These conditions lead to a simpler system, but such system describes no system of two 
coupled fields anymore. 
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We see that none of the conditions (0)"dZl) l ea ds to systems of coupled fields. How- 
ever, Eq. shows another possibility, which appears as follows: we investigate the 
quantity 

(Hw - H xx ) +4F^ X (78) 

with the aim of eliminating the square root in Eq. (|69D , restricting the set of parameters 
that defines H((j),x)- This possibility is particular, and can only be implemented in 
specific systems, where H(<fi, x) is explicitly known. For instance, in Ref. |23j one uses this 



reasoning to present analytical calculations that allow obtaining the discrete eigenvalues 
corresponding to BPS solutions of systems like the ones described by ([H]) or (0), in a 
specific region of the space of parameters. 

The above results show that although the BPS solutions are stable, in general it is 
not ease to know the corresponding eigenvalues explicitly. Evidently, we can always give 
up the analytical procedure to solve Schrodinger-like equations like Eq. fl70D numerically, 
but this is out of the scope of the present work. From the point of view of defects inside 
defects, the important information is that the BPS solutions are stable configurations, 
and so they can be considered stable defects to host non-BPS defects. The stability of 
the non-BPS defects will be investigated in the next subsection. 

The result that the general system does not give any system of two coupled fields 
for = H xx can be generalized in the following way. We introduce the new fields 
<p± = (0 ± x) an d use the equations of motion ([5]) and (|) to obtain 



dx 2 
dx 2 



H + H ++ + H_H_ + , (79) 
H + H + _ + EH . (80) 



Here we are using the notation H + = dH/d<p+, f/_ = dH/d<p-, and so forth. Derivatives 
of the function H are related by 

H<t> = "4 (H+ + H-), (81) 



V2 

H x = ^ (#+-#_), (82) 



and also 



1111 

H H = 2 H++ + 2 H ~ + 2 H+ ~ + 2 H + ' '~ 83 '' 

H xx = \h ++ + 1 -H^- 1 -H + _- 1 -H_ + , (84) 

H^ x = \h ++ - 1 -H^+ 1 -H_ + - 1 -H + _ , (85) 
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H x4> = i + (86) 

We use ( [S5l) and (^) to see that the condition H^ x = H x ^ now becomes H + _ = H_ + . 
Furthermore, from (|83|) and fl34|) we realize that if one further imposes the condition 
= H xx one gets = H_ + = 0. This result is general, and implies that H((j) + , 0_) 
can be written as the sum of two functions, one depending on + and the other on 0_. 
The Lagrangian density is then reduced to a sum of two Lagrangian densities, one for the 
field + and the other for 0_. The system decouples into two systems of a single field 
each one, and so it does not describe two coupled fields anymore. 



3.1 Stability of the topological solutions 

The pair of solutions with = and x given by Eq. ( PB| ) constitutes a pair of non- 
BPS solutions, that is, it does not obey the corresponding first-order equations. Then its 
stability should be investigated explicitly, and we do it now since it was not done in the 
former works 0, HI- Here, we have that U$ x vanishes at the classical values with = 0, and 
we are left with the following uncoupled Schrodinger-like equations, after appropriately 
rescaling the space coordinate, 



d 2 



+ 



dz 2 r\i 2 a 2 
( d 2 



- 4 + 2(2 + r) sech 2 zj r] n (z) 

,2 \ 



dz 2 ^ rfi 2 a 2 



4 + 6 sech z ( m ( z ) 








(87) 



which are well known modified Posch- Teller equations |24| whose eigenvalues are given by 

2 



r/j, 2 a 2 



rjj?a 2 



4 
4 



'2(2 + r) 



1 



1 

n H — 
2 



m + 



(89) 
(90) 



where n — 0, 1, .., < y2(2 + r) + 1/4—1/2 and m = 0,1. Since r > we find the following 
restriction on the solution ([33]) of the parity preserving model Hi 



< r 



< 1 . 



(91) 



However, since we assume that r ^ 1, in order to keep the fields coupled, we have 
r G (0,1) as the range of values that ensure stability of the non-BPS solution (|33|). 
Similar conclusions can be found for the other system, with Hi- 
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3.2 Stability of the periodic solutions 




We now discuss the stability of the periodic solutions (|54|) and (^). The general discussion 
on linear stability already presented applies equally well here, the difference in the present 
case is that the equations that appear for the fluctuations are now of the Lame type 

N(N + l)sn 2 (^)| f(z) = . (92) 

As we are going to show, some of these equations are solved by the Lame polynomials, 
since N can be identified with integer. But there are others for which N is not integer 
and so the exact solutions will be not completely known ||25| . Evidently, for each solution 
we have two equations of stability, the first describing fluctuation in (f) and the second in 
X- For the periodic solution (|54D the stability equations are 

-6A; 2 sn 2 (2)| rj(z) = , (93) 

1 + l)k 2 sn 2 (z)\((z)=0, (94) 
r r J 

where we omit indices on eigenf unctions and eigenvalues and also on k\ and ft2, for sim- 
plicity. The first stability equation is of the Lame type with N = 2, so that there are 
2N +1 = 5 solutions whose eigenvalues are 

w 2 




r 2 a 2 jji 2 b 2 (k) 



3, 3k 2 , 0, 1 + k 2 ± 2Vl -k 2 + k A . (95) 



Here, instability appears because of the last eigenvalue with the minus sign, which is 
always negative for any k since < k 2 < 1. 

The second stability equation is also of Lame type but now N = 2/r is not an integer, 
in general. The solutions for this case are not completely known but are expressible as 
series of elliptic functions which truncate to a polynomial when N becomes an integer. 
Also, the corresponding eigenvalues for h are not known in general, except for a few cases 
like N = 1/2 or 3/2. So, for special choices of the ratio r 

4 2 4/3 1 2/3 ... 



N 



1/2 13/2 2 3 



we are able to find explicit solutions and the corresponding eigenvalues. Let us analyse 
some of these possibilities. For r = 2, N = 1 the eigenvalues are: 



Aa 2 n 2 b 2 (k) ' 1 + k' 1 + k 
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which are non- negative except for the last one. For the case r = 1, N — 2 the eigenvalues 
are 



u 2 



., - = 3,3k 2 ,0,l + k 2 ±2Vl-k 2 + k 4 , (97) 

a 2 fi 2 b 2 (k) 

which are positive except for the last one with the minus sign, which is negative for any 
k G (0, 1). Let us also comment on one case where N in half-integer: we choose it to be 
N = 1/2, which corresponds to r = 4; in this case the eigenvalue is zero, exactly. 

A similar situation appears for the study of the stability of the second periodic solution 
(|55|), for which the stability equations are 

{& 2 + {1 + k2) (d^v + 6psn2 ^)} ^ = > ( 98 ) 
{& 2 + {1 + k2) (d^v + 1) - 2(2 + r ^ sn2 ^} rt*) = ■ (") 

Here we also find that N = 2 for the first stability equation, and in this case the eigenvalues 
are 

LO 2 



, - = 3, 3k 2 , 0, 1 + k 2 ± 2Vl - k 2 + A; 4 . (100) 
a z /j z rb z {k)k z 

Analogously to the previous case, the second stability equation for the solution (]55f) implies 

N = - l - + yi + S{2 + r) (101) 
so that in general we will have no closed solution except for special values of the ratio r 



r 


1 


4 8 


13 ... 


N 


2 


3 4 


5 ... 



Note that lower values of iV are not allowed since we imposed r > 0. However, if we 
choose r = 1, N = 2 we have the eigenvalues 



J 2 



a 2 ^ 2 b 2 (k) 



3, 3k 2 , 0, 1 + k 2 ± 2Vl - k 2 + k A , (102) 



and the periodic solution (|55|) is also unstable, as expected. 

We notice that the above stability equations recover the corresponding stability equa- 
tions of the former subsection in the limit k — > 1. Since sn 2 ^ — > tanh 2 z = 1— sech 2 ;?, we 
see that the operator 

d 2 

— + h k -k 2 psn 2 z, (103) 
dz l 
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changes to 



dz 2 



+ (hi — p) + p sech 2 2; . 



(104) 



However, there are subtleties in the limit k — *■ 1 for the eigenfunctions and eigenvalues of 
the corresponding periodic and topological solutions. For instance, it is only when k = 1 
that we have BPS solutions; for k 1 the periodic solutions will never solve the first-order 
equations. See Refs. [17|, [L8| for more detais on this issue in the case of a single scalar 
field. 



3.3 Energy and stability 

In this subsection we introduce another reasoning, that leads to another condition to be 
fulfilled by the classical solutions in order to make them stable against decaying into less 
energetic solutions. This reasoning is simple and can be interesting, mainly when one is 
unable to implement the standard investigation of stability analytically. The reasoning is 
based on the presence of stable BPS solutions, which solve the first-order equations, and 
so it does not work for periodic solutions. 

Let us consider the parity preserving model defined with Hi. In this case the non-BPS 
solution obtained with = has energy E = (4/3)/xr 3 / 2 a 3 , and connects the two vacuum 
states (0, a-s/r) and (0, —a^/r). However, these two vacuum states can be also connected 
by considering BPS solutions that make use of the vacuum (a, 0) or (— a, 0). Evidently, 
there are two degenerate possibilities of connecting the vacua (0, a-y/r) and (0, — ay/r): 
one uses (— a, 0) and the other (a, 0), choosing the left and right path, respectively. The 
energy of these two BPS solutions can be calculated easily - see Sec. |l]. The result is 
(4/3)/ira 3 . We compare this with the energy E = (4/3)/ir 3 / 2 a 3 of the non-BPS solution 
to see that it is only for r < 1 that the non-BPS defect does not decay into a pair of BPS 
defects. Interestingly, this is the same result we have already obtained using the standard 
investigation of classical or linear stability. 



4 High temperature effects 

Owing to the possibility of applications to cosmology, let us now study the effective poten- 
tial at finite temperature for the models introduced in this work. The effective potential 
allows introducing the conditions for symmetry restoration, to identify the symmetric 
phase, the phase in which the system supports no topological defects anymore. Here the 
investigations are done in the (3+1) dimensional space-time, and deal with constant and 
uniform field configurations. The one-loop effective potential and the corresponding finite 
temperature effects can be calculated according to the standard investigations |2(|. It can 
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be written as, keeping only the high temperature contributions ||, 

U T fa x) = Ufa X ) + + U xx ) • (105) 

The above expression for the effective potential shows that the thermal corrections add 
to mass terms in the following form 

m 2 H {T) = U^(0,0) + ^(U HH + U xx ^) , (106) 

m xx( T ) = U xx (0,0) + ^(U Hxx + U xxxx ) , (107) 

™l x (T) = U^(0,0) + ^(U^ x + U xxx ^ . (108) 

It is convenient to investigate the general model, defined with the potential given by 
Eq. (|i~3|). The effective potential at finite temperature in this case can be written as 

U T fa X ) = ^(A 2 + z/ 2 ) ( /) 4 + 2z/(A + yU )0 3 x + (A/i + 2^ + az/ + 2 / i 2 )0V 
+2fi(u + a)0 X 3 + \{f + a 2 ) X A + \nfo{T)<f + \ml x {T)x 2 
+ml x (T)<f>X + y A K,(0) + m 2 x (0)] + i(AV + a 2 6 4 ) , (109) 
where the mass parameters at finite and zero temperature are given by 





= m^(0) + ^-(3A 2 + 5z/ 2 + A yU + az/ + 2/i 2 ) , 


(110) 


-xx( T ) 


T 2 

= m 2 x (0) + — (3a 2 + A/i + 2z/ 2 + ( Tz/ + 5 y u 2 ) , 


(111) 




T 2 

= m 2 x (0) + — (A^ + (7// + 2H , 


(112) 


m xx(°) 


= -2(AV + afeV) , 
= -2(Aa 2 /i + crV) , 
= -2(Aa 2 z/ + crfc 2 //) . 


(113) 
(114) 
(115) 



In order to ensure that the potential mantains the Z 2 x Z 2 parity symmetry at finite 
temperature we impose, besides the restrictions already found at zero temperature, 

\v + afi + 2fiu = 0. (116) 
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We notice that this condition is fully satisfied by the choices ^ = a~ = 0orA = /z = 
corresponding to Hi or H 2 , respectively. This result indicates that the existing solutions of 
second-order equations at zero temperature may persist when the system is in equilibrium 
with a thermal bath. However, we have been unable to write the finite temperature 
effective potential U T ((p,x) via some generalized function H T ((j),x) in the form (f|) that 
appears at zero temperature. Although we have no explicit proof, we argue that the 
thermal effects destroy the Bogomol'nyi bound, together with the possibility of finding 
BPS defects at finite temperature. 

To define the masses properly and to introduce the critical temperatures we rotate 
the plane (0, x) to the plane (</>+, <p-) that diagonalizes the mass matrix. In this case we 
have m+(T) > and iriL(T) > as the conditions for symmetry restoration in each one 
of the two independet field directions <ft + and 0_, respectively. The critical temperatures 
are obtained in the limit where these masses vanish, and here they are given by 



c\2 . 12 



(miA 2 + m 2 Ai - 2m 3 A 3 ) ± v 7 ^ 



;il7) 



±! AiA 2 - A| 
where we have set 

A = (miA 2 - m 2 Ai) 2 + 4m l m 2 Al - 4m 3 (miA 2 A 3 + m 2 AiA 3 - m 3 AiA 2 ) . (118) 

The parameters m; and A 4 are derivatives of U(4>,x) evaluated at the point (0,0), and 
are given by 

mi = XJ H , m 2 = U xx , m 3 = , 

Ai = + U Hxx , A 2 = U xxxx + U^ xx , A 3 = U H(j)X + U xxx<p . (119) 

Let us consider for instance a model defined by the general potential U(4>, x) given by 
eq. (p~3| ) but with the restrictions b 2 = ra 2 , v = ra, and A = r/i, so that its vacuum states 
are (±a, 0) and (0, ±ay/r). The critical temperatures in this case are given by 



Y2a 2 r{r + l)[(r + 2)a 2 - (2r + l)/x 2 ] 
± ' " (5r + l)(2r 2 + r + 3)a 2 - (r + 5)(3r 2 + r + 2)^ 2 



c\2 



12a 2 r(r - l)^(r - l) 2 (a 4 + /x 4 ) + (6r 2 + 4r + 6)/x 2 (J 2 
(5r + l)(2r 2 + r + 3)a 2 - (r + 5)(3r 2 + r + 2)/z 2 



(120) 



Then, to find the critical temperatures for the parity preserving model Hi we just take 
a = on the above result so that we have 

12r 2 o 2 

ra 2 = o\l\ > ( m ) 

Sr z + r + 2 

iTlf = (122) 
r + 5 
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This result coincide with the one recently found in Ref. ||. 

It is also interesting to compute the critical temperatures for the model defined by H 2 , 
Eq. (|22|). They are also obtained from Eq. ( |120| ), but now imposing fi = 0. Here we get 



TO 2 = (123) 

(T c f = —P^ • (124) 

y 2r 2 + r + 3 V ' 

We notice that these results for H 2 are compatible with the ones obtained for Hi since 
they are connected via the identifications: ra 2 <-> a 2 and r <-> 1/r. 

The existence of two critical temperatures is due to the fact that the systems we are 
studing have two degrees of freedom and as a consequence their symmetries are restored 
separately for each one of the two independent field directions. 



5 Comments and conclusions 

In this paper we have dealt with systems of two coupled real scalar fields, searching for 
and investigating the corresponding classical or linear stability of the basic topological 
and periodic solutions when one of the two fields is set to zero. Our investigation is 
related to the possibility of nesting defects inside defects in 3 + 1 spacetime dimensions. 
However, in 3 + 1 dimensions renormalization restricts interactions to the fourth power 
in the polynomial potential, and so we have been mostly concerned with kinks of the 
hyperbolic tangent type. Such kinks have the property of vanishing at its own core or 
center, and so at the core of the host defect the system is reduced to a system of a single 
field, which is still able to generate the other defect, the defect to be nested inside the 
host defect. 

We have shown that the Z 2 x Z 2 parity symmetry is necessary to the formation of 
defects inside defects of this kind. As we have seen, the two basic defects that appears 
when one of the two fields vanishes can be summarized as follows: 



model 


symmetry 


topological 
solutions 


nontopological 
solutions 


H x 


Z 2 x Z 2 


(tanhx,0)/BPS 
(0, tanh x) / non-BPS 


(sn x, 0) 
(0, sn x) 


H 2 


Z 2 x Z 2 


(tanh x, 0) / non-BPS 
(0,tanhx)/BPS 


(sn x, 0) 
(0, sn x) 


E x or H 2 
(r = l) 


z 4 


(tanhx,0)/BPS 
(0,tanhx)/BPS 


(sn x, 0) 
(0, sn x) 
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We recall that in the last case, for r = 1 the two Zi x Zi models colapse to the same Z4 
model, but this model does not describe a system of two coupled fields anymore ||. 

We have also investigated the presence of periodic solutions when one of the two fields 
is set to zero. These solutions are direct extensions of the sphalerons found in |T7] in the 
case of a single real scalar field. The basic motivation for doing this is to enlarge the 
scope of the paper, since now one can mix topological and periodic defects, for instance 
allowing sphalerons to be nested inside domain walls. Furthermore, since the systems we 
have considered can be seen as real bosonic portions of supersymmetric models, it seems 
also interesting to investigate the behavior of fermions in such models. 

The high temperature results show the presence of two critical temperatures, signalling 
symmetry restoration in each one of the two independent field directions. Here the thermal 
effects seems to destroy the Bogomol'nyi bound, making the BPS defect to vanish at 
finite temperature. Within the standard cosmological evolution the temperature decreases 
until breaking the symmetry for the first field, allowing the formation of the host defect 
and later, after the breakdown of the symmetry corresponding to the second field, the 
second defect appears nested inside the first defect. As we can see, it is possible to 
built models in which topological defects can nest nontopological sphalerons, introducing 
instability inside the planar region on the wall, which may be of interest to applications to 
cosmology. Other interesting directions consider systems of coupled fields to investigate 
new aspects of inflation |^7|, |28|, ^UJ and the possibility of having topological defects 
ending on topological defects |3l], |32| . These and other related issues are presently under 
consideration. 
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